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Information geometric techniques and inductive inference methods hold great promise for solving 
computational problems of interest in classical and quantum physics, especially with regard to com- 
plexity characterization of dynamical systems in terms of their probabilistic description on curved 
statistical manifolds. In this article, we investigate the possibility of describing the macroscopic 
behavior of complex systems in terms of the underlying statistical structure of their microscopic de- 
grees of freedom by use of statistical inductive inference and information geometry. We review the 
Maximum Relative Entropy (MrE) formalism and the theoretical structure of the information ge- 
ometrodynamical approach to chaos (IGAC) on statistical manifolds Ms- Special focus is devoted to 
the description of the roles played by the sectional curvature ICmsi the Jacobi field intensity Jms 
and the information geometrodynamical entropy Sms (IGE). These quantities serve as powerful 
information geometric complexity measures of information-constrained dynamics associated with 
arbitrary chaotic and regular systems defined on A^5. Finally, the application of such information 
geometric techniques to several theoretical models are presented. 

PACS numbers: Probability Theory (02.50.Cw), Riemannian Geometry (02.40.Ky), Chaos (05. 45. -a), Com- 
plexity (89. 70. Eg), Entropy (89.70.Cf) 



I. INTRODUCTION 



It is commonly accepted that one of the major goals of physics is modeling and predicting natural phenomena using 
relevant information about the system of interest. Taking this statement seriously, it is reasonable to expect that 
the laws of physics should reflect the methods for manipulating information. Indeed, the less controversial opposite 
point of view may be considered where the laws of physics are used to manipulate information. This is exactly the 
point of view adopted in quantum information science where information is manipulated using the laws of quantum 
mechanics [H. In this work, we wish to explore an alternative viewpoint: perhaps the laws of physics are nothing 
but rules of inference . In this view the laws of physics are not laws of nature but are merely the rules we follow 
when processing the information that happens to be relevant to the physical problem under consideration. When the 
information available is sufficient to make unequivocal, unique assessments of truth we speak of making deductions: 
on the basis of this or that information we deduce that a certain proposition is true. In cases where we do not 
have statements that lead to unequivocal conclusions, we speak of using inductive reasoning, and the system for this 
reasoning is probability theory [3|. The word "induction" refers to the process of using limited information about 
a few special cases to draw conclusions about more general situations. Following this line of reasoning, interesting 
probabilistic approaches describing complex dynamical systems have already been investigated 0-01 ■ 

The study of the relationship between entropy and the complexity IT] of the trajectories of a dynamical system 
has always been an active field of research [8l-[lO|. Chaotic behavior is a particular case of complex behavior and it 
will be the subject of the present work. In this article we make use of the so-called Entropic Dynamics (ED) [Tl|. 
ED is a theoretical framework that arises from the combination of inductive inference (Maximum Entropy Methods 
(ME), [iMl]) and Information Geometry (IG) The most intriguing question being pursued in ED stems from 
the possibility of deriving dynamics from purely entropic arguments. This is clearly valuable in circumstances where 
microscopic dynamics may be too far removed from the phenomena of interest, such as in complex biological or 
ecological systems, or where it may be unknown or perhaps even nonexistent, as in economics. It has already been 
shown that entropic arguments do account for a substantial part of the formalism of quantum mechanics, a theory that 
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is presumably fundamental [l8l]. Perhaps the fundamental theories of physics are not so fundamental; they may just 
be consistent, objective ways of manipulating information. Following this line of thought, we extend the applicability 
of information geometric techniques and inductive inference methods to computational problems of interest in classical 
and quantum physics. In particular, we focus on the complexity characterization of dynamical systems in terms of 
their probabilistic description on curved statistical manifolds and identify relevant measures of chaoticity associated 
with such an information geometrodynamical approach to chaos (IGAC) 0-0, IT9l - [2^ . 

The layout of this article is as follows. In Section II, we review the Maximum relative Entropy formalism. In Section 
III, we present an introduction to the main features of the IGAC. In Section IV, we introduce the information-geometric 
indicators of chaos for our theoretical model: sectional (ICms) ^^.d Ricci (TZms) curvatures, Jacobi field intensity 
(Ja4s), IGE (Sms) s-n^d information geometric complexity (Cms)- In. Section V, we present five applications of our 
techniques to study the dynamical complexity of suitable statistical models. First, we characterize the chaotic behavior 
of an ED Gaussian model describing an arbitrary system of I degrees of freedom in absence of correlations. Second, 
we discuss the asymptotic temporal behavior of the information geometric complexity of the maximum probability 
trajectories of a 2/-dimensional Gaussian statistical manifold in presence of correlations between the macrovariables 
labeling the macrostates of the system f23'|. Third, an information-geometric analogue of the Zurek-Paz quantum 
chaos criterion of linear entropy growth i2Ql for a system of uncoupled, anisotropic (random frequency) and inverted 
Harmonic Oscillators is examined. Forth example, we present the information-geometric characterization of regular 
and chaotic quantum energy level statistics [l,!^- Our fifth and final example applies the IGAC to model and describe 
the scattering induced quantum entanglement of two micro-correlated, spinless, structureless, non-relativistic particles, 
where each particle is represented by a minimum uncertainty wave-packet (23 ) . Final remarks are presented in Section 
VI. 



II. ON THE MAXIMUM RELATIVE ENTROPY FORMALISM 

In 1957, Jaynes [25[ showed that maximizing statistical mechanic entropy for the purpose of revealing how gas 
molecules were distributed was equivalent to the maximizing of Shannon's information entropy [26j with statistical 
mechanical information. In traditional statistical mechanics, Boltzmann applied relevant information regarding the gas 
molecules in a closed system at equilibrium, whereas Gibbs used information relevant to an open system at equilibrium. 
However, as Jaynes noted, both solutions can be shown to be special cases of the maximum entropy principle. It is 
important to emphasize that this method of maximizing entropy is true for assigning probabilities regardless of the 
information specifics [l6| . As a consequence of this, the fact that Boltzmann and Gibbs used information relevant 
at equilibrium does not in general define entropy as only applying at equilibrium. Either of them could have instead 
included information that was time dependent as well. It is important to observe that Boltzmann included many 
assumptions into his solutions, such as assuming that all microstates were independent and had equal probability of 
occurring. The result of this was the interpretation that his probability assignments were for individual microstates. 
Gibbs on the other hand made no such assumptions. Therefore, his probability assignments described the system as 
a whole. This makes Gibbs method more general. In fact, it has been shown that Gibbs' entropy, not Boltzmann's, 
is numerically equivalent to thermodynamic entropy (27j . 

Jaynes' use of the method of maximum entropy for assigning probabilities is commonly known as MaxEnt '2?']. 
This method has evolved to a more general method, the method of Maximum (relative) Entropy (MrE) ilS-ilGi] which 
has the advantage of not only assigning probabilities but updating them (we use the term "updating" because as we 
gain valuable information, we "update" to new probabilities) when new information is given in the form of constraints 
on the family of allowed posteriors. This is similar in function to Bayes theorem, but using macroscopic information 
(expectation values) as opposed to microscopic information (data). We point out that one of the drawbacks of the 
MaxEnt method was the inability to include data. Although we do not include information in the form of data in 
this article, this problem has been solved by showing that the MrE method can be used when either data, constraint 



information or both is present [28|. For our purposes here, we will proceed to describe the simple case of how MrE is 
used to update from one probability to another given macroscopic constraint information only. Our first concern when 
using the MrE method to update from a prior to a posterior distribution is to define the space in which the search 
for the posterior will be conducted. Normally, when using the MrE method, the updating is based on three pieces of 
information: prior information (the prior) , the known or assumed relationship between the sets of microvariables and 
macrovariables (the model), and constraints on what the new model or posterior is allowed to be. Some good general 
examples of this are located in (29l.l30j. However, for our purposes in the present work we follow the simpler approach 
where our space of inquiry is limited to one dimension of one or several quantities, X £ X . We intend to maximize 
an appropriate entropy form, 

S[P, Poid] ^- j dX P{X) In (1) 
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subject to the appropriate constraints where P{X) = P{xi,X2, . . . ,xi), I is the dimensionahty of the microspace X 
and Poid (X) contains our prior information. We impose the usual normalization constraint, 

jdXP{X) = l, (2) 

and include additional information about X in the form of a constraint on the expected value of some function / (X), 

dXP{X)f{X)^{f{X))=F. (3) 



We proceed by maximizing ([T]) subject to constraints ([2]) and The purpose of maximizing the logarithmic relative 
entropy S[P, Poid] is to determine the value of P {X) that is closest to Pom {X) given the normalization and information 
constraints. Using the Lagrange multipliers formalism, we set the variation of 5'[P,Poid] with respect to P equal to 
zero, to obtain 



5{S[P,Po 



IdJ 



J dX P{X)-1 +(3 J dX P{X)f{X)-F ^ = 0. 



(4) 



The quantities a and /3 are Lagrange multipliers whose actual values are determined by the value of the constraints 
themselves. Substituting ([T]) into (jj]) we get (after some algebra) 



dX (in ^^^l +l-a-Pf{X)]sP (X) - 0, ^5P 

Therefore, the terms inside the curly brackets in Q must sum to zero, yielding 

P„cwW = PoidWe[-i+"+^^mi. 

In order to determine the Lagrange multipliers a and /3, we substitute our solution 
^ and ([3]), respectively. Substituting ^ into we obtain 



(5) 



(6) 



dXPoM (X)p/^W=e(i-"). 



P„ew(X) - Po\A {X) 



Then substituting ([7]) into ^ yields 



where 



The Lagrange multiplier (3 is determined by first substituting (jH]) into ([3]), 



- / dXe'^/^^^Poid {X) 



dX 



Now, P can be determined by rewriting pll|) as 



Poid {X) ■ 



dlnZ 



f{X)=F 



into the constraint equations 

(7) 

(8) 
(9) 

(10) 



dp 



= F. 



Finally, the selected posterior Pnow(-'^) has the expected canonical form 

e/3/(^) 

PncAX) = Pom{X)^—. 



(11) 



(12) 



When the expectation (j3|) is the first moment, the canonical form is clearly seen. However, if the second moment is 
used a normal distribution is produced, as in our example below. In conclusion, the ME method can be described 
as a general method to update from a prior distribution q{x) to a posterior distribution p{x) when new information 
becomes available. It is perhaps worth emphasizing that in this approach, entropy is a tool for reaso ning which 
requires no interpretation in terms of heat, multiplicities, disorder, uncertainty, or amount of information |16j. 

In what follows we will use this inference methodology together with information geometric methods to characterize 
the temporal complexity (chaoticity, dynamical stochasticity) of suitable statistical models of physical relevance. 
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III. ON THE INFORMATION GEOMETRY OF CHAOS 



IGAC is an application of ED to complex systems of arbitrary nature. It is the information-geometric analogue of 
conventional geometrodynamical approaches ^, [H, [3l| - [33| where the classical configuration space Te is replaced by a 
statistical manifold A4s i-^s is defined below in (dH])) with the additional possibility of considering chaotic dynamics 
arising from non conformally fiat metrics (the Jacobi metric is always conformally fiat, instead). It is an information- 
geometric extension of the Jacobi geometrodynamics (the geometrization of a Hamiltonian system by transforming it to 
a geodesic flow (33|). The reformulation of dynamics in terms of a geodesic problem allows for the application of a wide 
range of well-known geometrical techniques to the investigation of the solution space and properties of the equations 
of motion. The power of the Jacobi reformulation rests in the fact that all dynamical information is collected into a 
single geometric object (namely, the manifold on which geodesic flow is induced) in which all the available manifest 
symmetries of the system are retained. For example, the integrability of a system is connected with the existence 
of Killing vectors and tensors on this manifold. The sensitive dependence of trajectories on initial conditions, which 



is a key ingredient of chaos, can be investigated from the equation of geodesic deviation. In the Riemannian [31 1 
and Finslerian (s^ (a Finsler metric is obtained from a Riemannian metric by relaxing the requirement that the 
metric be quadratic on each tangent space) geometrodynamical approach to chaos in classical Hamiltonian systems, 
an active field of research concerns the possibility of finding a rigorous relation among the sectional curvature, the 
Lyapunov exponents, and the Kolmogorov-Sinai dynamical entropy (i.e., the sum of positive Lyapunov exponents) 
[35| . Using information-geometric methods, we have investigated in some detail the aforementioned open research 
problem [IB, d-il . 

IGAC arises as a theoretical framework to study chaos in informational geodesic fiows describing physical, biological 
or chemical systems 0-3 US] • The main goal of an ED model is that of inferring "macroscopic predictions" in the 
absence of detailed knowledge of the microscopic nature of the arbitrary complex systems being considered. More 
explicitly, by "macroscopic prediction" we mean knowledge of the statistical parameters (expectation values) of the 
probability distribution function that best refiects what is known about the system. This is an important conceptual 
point. The probability distribution refiects the system in general, not the microstates. We then select the relevant 
information about the system. In other words, we have to select the macrospace of the system. It is worth mentioning 
that the coexistence of macroscopic and microscopic dynamics for a given physical (biological, chemical) system from 
a dynamical and statistical point of view has always been a very important subject of investigation [36j . For example, 
in certain fiuid systems showing Rayleigh-Benard convection [37| . the macroscopic chaotic behavior (macroscopic 
chaos) is a manifestation of the underlying molecular interaction of a very large collection of molecules (microscopic or 
molecular chaos). Moreover, macroscopic chaos arising from an underlying microscopic chaotic molecular behavior has 
been observed in chemical reactions. In order to study the underlying dynamics in rate equations of chemical reactions, 
a mesoscopic description has been adopted which is given by a set of transition probabilities among chemicals. In 
such a description, the underlying dynamics of macroscopic motion is that of stochastic processes and the evolution 
of the probability distribution associated with each chemical is investigated [33i] . 

In what follows, we schematically outline the main features underlying the construction of an arbitrary form of 
entropic dynamics. First, the microstates of the system under investigation must be defined. For the sake of simplicity, 
we assume the system is characterized by an ^-dimensional microspace X with microstates X = {xi, . . . ,xi). The 
quantities Xk with k = l,2,...,l denote the degrees of freedom of the microstates of the system. We assume 
that each degree of freedom Xk is subject to n-information constraints. The n-information constraints may (though 
not necessarily) correspond to the moments of P {X\Q). For example, if we consider a system in which the relevant 
information constrains are given by the first two moments of P {X\Q), namely the expectation values /ij, and variances 



{Xk) 



Mfe and (^(^{xk - (a^fc))^)) ^ = CTfe, (13) 



then the resulting posterior probability distribution function is a Gaussian distribution pk (xfej/i^,, ak)- For general 
purposes, it is convenient to define an n/-dimensional macroscopic vector Q with statistical coordinates (where 

k = 1,2, . . . ,1 labels the microstates and m = 1,2, ... ,n enumerates the information constraints). The set Q defines 
the nZ-dimensional space of macrostates of the system, namely the parameter space Iq . In the Z-dimensional Gaussian 
case, the vector Q is n/-dimensional and has coordinates {^i, . . . , ^i;ai, . . . ,ai). 

At this point we make a working hypothesis that the microstates are statistically independent. Then, in addi- 
tion to information constraints, each distribution pk {^k\"d^k^^^ of each degree of freedom Xk must satisfy the usual 
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normalization conditions 



dxkPk [xkl^i""^) = 1- (14) 



Once the microstates have been defined and the relevant (linear or nonlinear) information constraints selected, we are 
left with a set of Z-dimensional vector probability distributions 

p{x\e) = f[pk{x,\4"^^) (15) 

fe=i 

(assuming statistical independence, identically distributed microvariables and assuming as an additional working 
hypothesis that the prior probability distribution is uniform) encoding the relevant available information about the 
system. 

The statistical manifold Ms, 

Ms - |p(^|e) = flpk (16) 

is defined as the set of probabilities {P {X\Q)} described above with X G M}, 6 G Xe where Xq is comprised of the n 
parameter sub-spaces corresponding to each of the n-information constraint according to Xq — x'™^ ['^k"l ■ Each 

m=l ^ ^ 

point of the geodesic on an n?-dimensional statistical manifold Ms represents a macrostate parametrized by the 
macroscopic dynamical variables defining the macrostate of the system. Furthermore, each macrostate is in 

one-to-one correspondence with the probability distribution P {X\Q) representing the maximally probable description 
of the system being considered. Thus, the set of macrostates forms the parameter space le while the set of probability 
distributions forms the statistical manifold Ms- Referring to our Gaussian example above, Xq is the direct product 
of the parameter sub-spaces X^ (corresponding to the first moment n = 1 of the Gaussian, i.e., the expectation /i) 
and Xcr (corresponding to the second moment n = 2 of the Gaussian, i.e., the variance a), where 1^ = (— oo,+oo)^ 
and Icr = (0, -l-oo)^. 

A measure of distinguishability among macrostates is obtained by assigning a probability distribution P {X\0) 3 
A^5 to each macrostate 9. Assignment of a probability distribution to each state endows Ms with a metric structure. 
Specifically, the Fisher- Rao information metric gab (©) flTf , 

gabie)= J dXP{X\e)dalnP{X\Q)db\nP{X\Q); a,b^l,...,nl andda^ (17) 

defines a measure of dis ting uishability among macrostates on A^s- 

It is known from IG [T3| that there is a one-to-one relation between elements of the statistical manifold Ms and 
the parameter space Xq. Specifically, the statistical manifold A^s is homeomorphic to the parameter space Xq. This 
implies the existence of a continuous, bijective map Hm^.Tq, 

hM.,Xe ■■Ms3 P{X\e)^e£XQ, (18) 

where hj^ (0) = P {X\Q). The inverse image hj^ is the so-called homeomorphism map. It is worth pointing 
out that the possible chaotic behavior of the set of macrostates Q is strictly related to the selected relevant information 
about the set of microstates X of the system. In other words, the assumed Gaussian characterization of the degrees 
of freedom {xk} of each microstate of the system has deep consequences on the macroscopic behavior of the system 
itself. More generally, within our theoretical construct, "the macroscopic behavior of an arbitrary complex system is a 
consequence of the underlying statistical structure of the microscopic degrees of freedom of the system being considered". 

It should be noted that coupled constraints would lead to a "generalized" product rule in ([T5|) and to a metric 
tensor (|17p with non-trivial off-diagonal elements (covariance terms). In presence of correlated degrees of freedom 
{xj}, the "generalized" product rule becomes 

n n 
-Ptot ( ) = \IP3 i^J ) Ptot ixi,...,Xn)^Y[P, {Xj ) , (19) 

i=i i=i 
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where 

Pni ■^n 1^1 ; • ■ • ; •^n— 1 )Pn-l{ 2) ■ ■ ■ P2{x2\xi) Pi{xi) . (20) 

For instance, correlations in the degrees of freedom may be introduced in terms of the following information-constraints, 

= fj i^i^ ■ • ■ ' ^3-1) ' = 2, . . . , n. (21) 

In such a case, we obtain 

Ptot (Xl,- ■ • ,a;„) = 6 {Xn - /„ (Xi, . . . ,a:„_i)) S (Xn^l - fn-l {Xl, ■ ■ .,Xn-2)) ' ' ' S {X2 ~ f2 (xi)) Pi (xi) , (22) 

where the j-th probability distribution Pj (xj) is given by 



P] (Xj) ^ J ■ ■ ■ J dxi ■ ■ ■ dXj-ldXj + l ■ ■ ■ dxnPtot (xi,. . . ,Xn) ■ 



(23) 



Correlations between the microscopic degrees of freedom of the system {xj} (micro-correlations) are conventionally 
introduced by means of the correlation coefficients r. 



(micro) 

r- ■ 



T (xi , Xj ) 



(micro) 



dcf {XiXj) {Xi) {Xj) 



{xi - {xi)Y 



(24) 



(micro) 



€ (—1, 1) and i,j = 1, . . . ,n. For the 2n-dimensional Gaussian statistical model in presence of micro- 



with r, 

correlations, the system is described by the following probability distribution P{X\Q), 



P{X\Q) 



[(27r)"detC(e)]- 



■ exp 



-- (X - M)* • C-^ (6) • (X - M) 



vrcr?) ^ exp 



2a| 



, (25) 



where X = (xi, . . . , x„), M — (/ij^, . . . , /i„) and C (O) is the {2n x 2n)-dimensional (non-singular) covariance matrix. 

Once AAs and Xq are defined, the ED formalism provides the tools to explore dynamics driven on AAs by 
entropic arguments. Specifically, given a known initial macrostate (probability distribution), and that the 

system evolves to a final known macrostate the possible trajectories of the system are examined in the ED 

approach using ME methods. We emphasize ED can be derived from a standard principle of least action (Maupertuis- 
Euler-Lagrange-Jacobi-type) .IL .38, ] . The main differences are that the dynamics being considered here, namely 
Entropic Dynamics, is defined on a space of probability distributions AAs, not on an ordinary linear space V and the 
standard coordinates qa of the system are replaced by statistical macrovariables i?". The geodesic equations for the 
macrovariables of the Gaussian ED model are given by nonlinear second order coupled ordinary differential equations 



rfr^ ^'^ dr dr 



(26) 



The geodesic equations in (pS)) describe reversible dynamics whose solution is the trajectory between an initial 0('nit'ai) 
and a final macrostate The trajectory can be traversed equally well in both directions. A geodesic on a 

curved statistical manifold Ms represents the maximum probability path a complex dynamical system explores in 
its evolution between initial and final macrostates Q ('initial) ^^^^ q (final) ^ respectively. We point out that this path is 
obtained by use of a Principle of Probable Inference, namely the Maximum relative Entropy Method. Our theoretical 
formalism allows us to analyze important physics problems through statistical inference and information geometric 
techniques. 



IV. INFORMATION GEOMETRIC MEASURES OF TEMPORAL COMPLEXITY 

In this Section, mainly following one of our previous works j2^ , we introduce the relevant indicators of chaoticity 
(temporal complexity, dynamical stochasticity; in general there is no one-to-one relation between chaos and complexity, 
chaos may imply complexity but necessarily vice versa) within our theoretical formalism. Specifically, once the 
Fisher- Rao information metric is given, we apply standard methods of Riemannian differential geometry to study the 
information-geometric structure of the manifold A4s underlying the entropic dynamics. Connection coefficients F^^, 
Ricci tensor TZab, Riemannian curvature tensor TZabcd, sectional curvatures ICms [39j . scalar curvature TZms: Weyl 
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anisotropy tensor Watcd (the anisotropy of the manifold underlying the system dynamics plays a significant role in the 
mechanism of instability), the Jacobi vector field intensity Jms^ the information-geometric analogue of the Lyapunov 



exponent A j and the information-geometric entropy (IGE) Sm, can be calculated in the standard manner [22ll3lLl32 

In order to characterize the chaotic behavior of complex entropic dynamical systems, we are primarily concerned 
with the signs of the scalar and sectional curvatures /Cwig of the asymptotic behavior of Jacobi fields J° on 
M.s^ the existence a non-vanishing Weyl anisotropy tensor Wafccrf and the asymptotic behavior of the IGE Sms- It is 
crucial to observe that true chaos is identified by the occurrence of two features |32| : 1) strong dependence on initial 
conditions and exponential divergence of the Jacobi vector field intensity, i.e. stretching of dynamical trajectories; 2) 
compactness of the configuration space manifold, i.e. folding of dynamical trajectories. 

A. Sectional and Scalar Curvatures 

Once the Fisher-Rao information metric gab is given, we use standard differential geometry methods applied to 
the space of probability distributions to characterize the geometric properties of A^s. Recall that the Ricci scalar 
curvature is given by 

■R = g"''"Rab, (27) 
where g°'^gbc = <^°c so that g°-^ — {gab)~^ ■ The Ricci tensor TZab is given by 

TZab = d^V:, dbV:^ + K,Ti, - TtTt,. (28) 
The Christoffel symbols rj^j, appearing in the Ricci tensor are defined in the standard manner as 

Kb = ^f"^ {dagdb + dbgad - ddgab) ■ (29) 
It can be shown that the Ricci scalar curvature can be written as (39j 

TZm. = TZabcdg^'g''' = 7^°„ - ^ /C (e„ e,) . (30) 

The scalar curvature is the sum of all sectional curvatures /C (e^, e,) of planes spanned by pairs of orthonormal basis 
elements {ca — 9^«(p)} of the tangent space TpMs with p e Ms [3a|, 

irl \ nabcdU°-V^u''v'^ i -^u-/- /qiN 

IC {u, V) = . a b c d '^ h , V ^h-' with I ^ J, (31) 

[gadgbc - gacgbd) u^v^u^'v" 

where {ea,h'>) = S\. Not ice that the sectional curvatures completely determine the curvature tensor. 

The negativity of the Ricci scalar IZms implies the existence of expanding directions in the configuration space 
manifold Ais- Indeed, since TZms is the sum of all sectional curvatures of planes spanned by pairs of orthonormal 
basis elements {cq = d^-^}, the negativity of the Ricci scalar is only a sufficient (not necessary) condition for local 
instability of geodesic flow. For this reason, the negativity of the scalar provides a strong criterion of local instability. 
Scenarios may arise where negative sectional curvatures are present, but the positive ones could prevail in the smn so 
that the Ricci scalar is non-negative despite the instability in the flow in those directions. Consequently, the signs of 
K-Ms '^^ crucial signiflcance for the proper characterization of chaos. 

B. Killing Vectors 

Yet another useful way to understand the anisotropy of the is the following. It is known that in N dimensions, 
there are at most '^'^'^^^'^ independent Killing vectors (directions of symmetry of the manifold). Since A^s is not a 
pseudosphere, the information metric tensor does not admit the maximum number of Killing vectors Ka defined as 

Cxgab = VaKb + VbKa = 0, (32) 

where 2?a, defined as 

VaKb = daKb - n^K,, (33) 
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is the covariant derivative operator with respect to the connection F defined in (l29l) . The Lie derivative CKQab of the 
tensor field gab along a given direction K measures the intrinsic variation of the field along that direction (that is, the 
metric tensor is Lie transported along the Killing vector) [iOl- Locally, a maximally symmetric space of Euclidean 
signature is either a plane, a sphere, or a hyperboloid, depending on the sign of TZ. In our case, none of these scenarios 
occur. As will be seen in what follows, this fact has a significant impact on the integration of the geodesic deviation 
equation on Ais- At this juncture, we emphasize it is known that the anisotropy of the manifold underlying system 
dynamics plays a crucial role in the mechanism of instability. In particular, fiuctuating sectional curvatures require 
also that the manifold be anisotropic. However, the connection between curvature variations along geodesies and 
anisotropy is far from clear . Krylov was the first to emphasize [4lj the use of 7?, < as an instability criterion in 
the context of an A'^-body system (a gas) interacting via Van der Waals forces, with the ultimate hope to understand 
the relaxation process in a gas. However, Krylov neglected the problem of compactness of the configuration space 
manifold which is important for making inferences about exponential mixing of geodesic fiows (42j . Compactness 
[43I lil j is required in order to discard trivial exponential growths due to the unboundedness of the "volume" available 
to the dynamical system. In other words, the folding is necessary to have a dynamics that is actually able to mix 
the trajectories, making practically impossible, after a finite interval of time, to discriminate between trajectories 
which were very nearby each other at the initial time. When the space is not compact, even in presence of strong 
dependence on initial conditions, it could be possible in some instances (though not always), to distinguish among 
different trajectories originating within a small distance and then evolved subject to exponential instability. As a 
final remark, we note that since homeomorphisms preserve compactness and since AAs and le are homeomorphic, it 
is sufficient to show that the parameter space Tq is compact in order to ensure that the corresponding submanifold 
of A^s is itself compact. 



C. Jacobi Fields 



A powerful mathematical tool to investigate the stability or instability of a geodesic flow is the Jacobi-Levi-Civita 

The JLC-equation covariantly describes how nearby geodesies 



equation (JLC equation) for geodesic spread [31 



locally scatter and relates the stability or instability of a geodesic flow with curvature properties of the ambient 
manifold. For the sake of clarity, consider the behavior of a family of neighboring geodesies {'&'ms ^)}fGR+« 
a statistical manifold M.s with dimaA^s — N. The geodesies iJ^U (''';?) are solutions of equation ((26)) . The relative 
geodesic spread on a non- maximally symmetric curved manifold as A^^ is characterized by the Jacobi-Levi-Civita 
(JLC) equation [H IH 



2?'^" dd' d^" _ 



(34) 



where the covariant derivative in ([M)) is defined as [4 



, ddf dr , 



„ d^"^ dr 
^"^"^ dr dr 



and the Jacobi vector field J" is given by [47 



(35) 



(36) 



Equation forms a system of N coupled ordinary differential equations linear in the components of the devia- 
tion vector field ([36)) but nonlinear in derivatives of the metric ([TT)) . It describes the linearized geodesic fiow: the 
linearization ignores the relative velocity of the geodesies. When the geodesies are neighboring but their relative 
velocity is arbitrary, the corresponding geodesic deviation equation is the so-called generalized Jacobi equation ps} . 
The nonlinearity is due to the existence of velocity-dependent terms in the system. Neighboring geodesies accelerate 
relative to each other with a rate directly measured by the curvature tensor TZatcd- The non-trivial integration of ()34p 
leads to the following expression of the Jacobi vector field intensity. 



rw'^'igabrj'y 



(37) 
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For applications of the asymptotic temporal behavior of Jms (''") ^ ^ reliable indicator of chaoticity, we refer to 
our previous articles in references Sllllil. The geodesic spread GO. M.S is described by means of an exponentially 
divergent Jacobi vector field intensity J"^, , a classical feature of chaos. In our approach, inspired by the work presented 
in [49[, the quantity A^^s defined as 



1 I \JMs{r)\'^ 
— In 



AJms ('^) 



dT 



dJMg (r) 




dT 


T = 



(38) 



characterizes the exponential growth rate of average statistical volumes (see ([55]) ) in A^^. This suggests that Xms 
may play the same role as the standard Lyapunov exponent (soj . Lyapunov exponents are asymptotic quantities since 
they are defined in the limit as time approaches infinity. 

D. Weyl Projective Curvature Tensor 

The Weyl Projective curvature tensor (or the anisotropy tensor) Wabcd is defined as 



hcd 



hcd 



N[N -I) 



{9bd9ac - gbc9ad) , 



(39) 



where N = nl is the dimension of the curved manifold. In (j39|) . the quantity TZatcd is the Riemann curvature tensor 
defined in the usual manner by 



bed — (^c^ bd ~ "(jr be + ^ /c^ bd~ ^ fd^ be 



(40) 



Considerations regarding the negativity of the Ricci curvature as a strong criterion of dynamical instability and the 
necessity of compactness of A^s in "true" chaotic dynamical systems requires additional investigation. 



E. Information Geometric Complexity 



Once the distances among probability distributions have been assigned, a natural next step is to obtain measures for 
extended regions in the space of distributions. Consider an TV-dimensional volume of the statistical manifold A^s of 
distributions P {X\S) labelled by parameters "i?" with a — 1, . . . ,N. The parameters i?" are coordinates for the point 
P and in these coordinates it may not be obvious how to write an expression for a volume element dVM^- However, 
within a sufficiently small region (volume element) any curved space looks fiat. Curved spaces are "locally flat". The 
idea then is rather simple: within that very small region, we should use Cartesian coordinates and the metric takes a 
very simple form, namely the identity matrix Sab- In locally Cartesian coordinates the volume element is simply 
given by the product 



which, in terms of the old coordinates is 

dVM^ = 



dx 



de 



dd^dd^ ■ ■ ■ d§^ = 



dx 



de 



(41) 



(42) 



ae 



of the transformation that takes the metric gab into its 



The problem at hand is the calculation of the Jacobian 
Euclidean form Sab- 

Let the new coordinates be defined by x° = S° . . . , . A small change dQ corresponds to a small change 
dx, 



dx"' = X^'dr where X^' "^^^ 



or 



and the Jacobian is given by the determinant of the matrix X^ , 



dx 




dct (x:') 


de 







(43) 



(44) 



10 



The distance between two neighboring points is the same whether we compute it in terms of the old or the new 
coordinates, 



Therefore, the relation between the old and the new metric is 



gab — Sa'b'X^ . 



Taking the determinant of p6)) . we obtain 



g =^ det (gab) 



det 



and therefore 



det [x:') I ^ V?- 



(45) 
(46) 

(47) 
(48) 



Finally, we have succeeded in expressing the volume element totally in terms of the coordinates Q and the known 
metric ffafc (6), 



The volume of any extended region on the manifold is given by 



Vm. = J dVM. - j Vad^'Q. 



(49) 



(50) 



Observe that y^d^Q is a scalar quantity and therefore is invariant under general coordinate transformations 0—^0', 
preserving orientation. The square root of the metric tensor transforms according to 



VgW) 



de 



and the flat infinitesimal volume element d Q transforms as 



d^e 



de 



56' 



Thus, from ((5T|l and ((52|) we obtain 



(51) 



(52) 



(53) 



Equation (j53p implies that the infinitesimal statistical volume element is invariant under general coordinate transfor- 
mations that preserve orientation, that is with positive Jacobian. 
The volume of an extended region AVm, {t; ?) of A^s is defined by 



e(r;?) 



(54) 



e(0;f) 



where 8 (r; s*) are solutions of the geodesic equations (|26p and — (ij^, . . . , <r^) is the quantity parameterizing 
the family of geodesies {i^^U^ (''"' ?)}^gR+N^- quantity that encodes relevant information about the stability of 
neighboring volume elements is the average volume (information geometric complexity) Cm^ (t; ?) defined as [l^ 



Cm. [r; r) {AVm. (r; ?)), - ^ J AVm. (r'; ?) dr' 





(55) 



We will call Cm. (t't?) the information geometric complexity of the maximally probable trajectories 0(t;^). The 
IGC Cm. (t'i ?) in (j55p represents the temporal average of the iV-fold integral over maximum probability trajectories 
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(geodesies) and serves as a measure of the number of the accessible macrostates in the configuration (statistical) 
manifold after a finite temporal increment r. In other words, {''']?) can be interpreted as the temporal 
evolution of the system's uncertainty volume Cm^ (0; For instance, Cm^ (0; ?) ni^y be a spherical volume of initial 
points whose center is a given point on an attractor and whose surface consists of configuration points from nearby 
trajectories. An attractor is a subset of the manifold Ms toward which almost all sufficiently close trajectories 
converge asymptotically, covering it densely as time goes on. Strange attractors are called chaotic attractors. Chaotic 
attractors have at least one finite positive Lyapunov exponent 52]. As the center of Cx, (0;<f) and its surface points 
evolve in time, the spherical volume becomes an ellipsoid with principal axes in the directions of contraction and 
expansion. The average rates of expansion and contraction along the principal axes are the Lyapunov exponents [53] . 



F. Information Geometric Entropy 



Finally, the asymptotic regime of diffusive evolution describing the possible exponential increase of average volume 
elements on provides another useful indicator of dynamical chaoticity. The exponential instability characteristic 
of chaos forces the system to rapidly explore large areas (volumes) of A^^- We remark that the exponential instability 
does not necessarily imply that trajectories explore large areas (volumes) of A^^- The folding mechanism is responsible 
for a dense exploration of areas (volumes) that may be small in principle. It is interesting to note that this asymptotic 
behavior appears also in the conventional description of quantum chaos where the von Neumann entropy increases 
linearly at a rate determined by the Lyapunov exponents. The linear increase of entropy as a quantum chaos criterion 
was introduced by Zurek and Paz (s^ l . In our information-geometric approach a relevant quantity that may be useful 
to study the degree of instability characterizing ED models is the information geometrodynamical entropy (IGE) 
defined as i-i, Ili42l 



Sm. (t; ?) = Sm. [e (r; ?) , 6 (0; ?)] = ^lim In | ^ ^ dr' J ^d^'S 
where g — |det {gab)\- More synthetically. 



e(T;?) 



e(o;<r) 



(56) 



dof 



(57) 



with (t;?) being the information geometric complexity defined in (|55|) . The IGE is intended to capture the 
temporal complexity (chaoticity) of ED models on curved statistical manifolds AJ^ by considering the asymptotic 
temporal behaviors of the average statistical volumes occupied by the evolving macrovariables labelling points on A^^. 



V. APPLICATIONS 



In this Section we present five applications of the IGAC. First, we study the chaotic behavior of an ED Gaussian 
model describing an arbitrary system of / uncorrelated degrees of freedom and show that the hyperbolicity of the 
non- maximally symmetric 2Z-dimensional statistical manifold Ms underlying such an ED Gaussian model leads to 
linear IGE growth and to exponential divergence of the Jacobi vector field intensity [l^]. Second, we study the 
asymptotic behavior of the dynamical complexity of the maximum probability trajectories on Gaussian statistical 
manifolds in presence of correlation-like terms between the macrovariables labeling the macrostates of the system 
under investigation. In presence of correlation-like terms, we observe a power law decay of the information geometric 
complexity at a rate determined by the correlation coefficient [23j . We also present an information-geometric analogue 
of the Zurek-Paz quantum chaos criterion of linear entropy growth [s^ . This analogy is presented by studying 
the information geometrodynamics of an ensemble of macroscopic, random frequency, inverted harmonic oscillators 
[1, HqI- Next, we apply the IGAC to study the entropic dynamics on curved statistical manifolds induced by classical 
probability distributions in common use in the study of regular and chaotic quantum energy level statistics. In doing 
so, we suggest an information-geometric characterization of regular and chaotic quantum energy level statistics [3, [I]- 
Finally, we apply the IGAC to characterize the quantum entanglement produced by a head-on collision between two 
Gaussian wave packets interacting via a scattering process i55,] . 
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A. Gaussian Statistical Model in Absence of Correlations 



As a first example, we apply the IGAC to study the dynamics of a system with I degrees of freedom, each 
one described by two pieces of relevant information, its mean expected value and its variance (Gaussian statisti- 
cal macrostates). The line element ds^ = gab (B) dd°'dd^ {a,h = 1, . . . , 21) on M.s is defined by [g^l 



k=l 



=2^[ ^^t'k + -^dai ) . (58) 



This leads to consider an ED model on a non-maximally symmetric 2Z-dimensional statistical manifold AAs- Manifold 
Ais possesses a constant negative Ricci curvature that is proportional to the number of degrees of freedom of the 
system, Rm, = The system explores statistical volume elements on Ais at an exponential rate, while the informa- 
tion geometrodynamical entropy Sm^ increases linearly in time (statistical evolution parameter) and is proportional 
to the number of degrees of freedom of the system, IXt. The parameter A characterizes the family of 

probability distributions on Alg. For the case being considered here, A does indeed play the role of the standard 
Lyapunov exponent fsoj . 

Recall that the finite Lyapunov exponent in the direction v € R^' of a trajectory (r, 9o) satisfying the differential 
equation Q = A{t)Q with 6 G M?'' and initial condition Q (0, 9o) = 9o is defined as |56j] 



A (w) '= lim In 



yJ{Xv,Xv) 



(59) 



The brackets (•, •) in (l59l) denote the standard scalar product in M^' and X = X {t]Q (r, 9o)) is the asymptotically 
regular fundamental matrix of the differential equation Q — A{t)Q = [s^l- For instance, in the 2/-dimensional 
Gaussian statistical model considered here, the set of differential equations to consider is 

^^-^(T)e(r)=0, (60) 

where (r) is the 2Z-dimensional vector Q (r) = {^^ (t) , . . . , /i; (r) ; ui (r) ^ . . . ,ai (r)) whose components are so- 
lutions of geodesic equations describing the evolution of the macrostates of the system [l^l- In the asymptotic 

limit, the 21 x 21 matrix A (t) can be approximated by a diagonal matrix with constant coefficients, A (r) w 
diag (0, Ai, 0, A2, . . . , 0, A;). A straightforward calculation leads to an asymptotically regular 21 x 21 fundamental 
matrix 

X (t) ^ diag (ciT, C2 exp (Air) , C3T, C4 exp (Asr) , . . . , c/_iT, c; exp (A/t)) , (61) 
with Ci e M, Vi = 1, . . . , L Therefore, equation ((59|) leads to the following interesting result, 

Amax (f ) = max { Ai , . . . , A, } , Vw e . (62) 

K+ 

Thus, the quantities Afc with k = 1,. . .,1 are indeed Lyapunov exponents. In this case, for the sake of simplicity, 
we have assumed Xi = Xj Wi, j — 1, . . . ,1. The asymptotic linear information-geometrodynamical entropy growth of 
(''') may be considered the information-geometric analogue of the von Neumann entropy growth introduced by 
Zurek-Paz, a quantum feature of chaos. 

The geodesies on Ai^ are hyperbolic trajectories. Using the Jacobi-Levi-Civita (JLC) equation for geodesic spread, 
we show that the Jacobi vector field intensity J^Ms diverges exponentially and is proportional to the number of degrees 
of freedom of the system, J^Ms /exp(Ar). The exponential divergence of the Jacobi vector field intensity J^Ms 

is a classical feature of chaos. Therefore, we conclude that [22] 

TZm, = —I, Jm, ^ exp (At) , Sm, IXt. (63) 

By virtue of we observe that TZmsj ^Ms ^^'^ ^Ms ^'''^ proportional to the number of Gaussian-distributed 
microstates of the system. This proportionality, even though proven in a very special case, suggest there may be a 
substantial link among these information-geometric indicators of chaoticity. 
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B. Gaussian Statistical Model in presence of Macro-correlations 



As a second example, we apply the IGAC to study the information constrained dynamics of a system with I 
degrees of freedom, each one described by two "correlated" pieces of relevant information, its mean expected value 
and its variance (Gaussian statistical macrostates). The line element ds^ = gab (©) dd°'dd^ (a, & = 1, ... , 21) on the 
2Z-dimensional Gaussian statistical manifold A^s in the presence of non trivial off-diagonal terms is given by 




We consider positive coefficients rj £ (0, 1), Vj — I, . . . ,1. From (|64|) . it can be shown that the Ricci scalar curvature 
Ti-Ms (''ii • ■ • I n) of such a 2?-dimensional manifold is given by 

k=l 

Notice that in the limit of vanishing correlation strengths {ru}, Ti-Ms = ~^ shown in [12]. Applying the IGAC 
formalism, we are able to compute the asymptotic temporal behavior of the dynamical complexity of geodesic trajec- 
tories for the correlated 2/-dimensional Gaussian statistical model. The technical details that will be omitted in what 
follows may be found in [2^. It turns out that [2^ 



21 



Sms (T;{Afc},{rfc}) 



El- 

k=l 



Ai (rfc) + 



A2 (rfe, Afc) 



(65) 



where 



1 + ^1 + Arl rk Afc 2 

The quantity S (r^, , Afc , a-t ) is a strictly positive function of its arguments. For = Vfc and s = 1,...,/, the 
information geometric entropy ^ A, r) becomes 



In 



Ai (r) + 



A2 (r,A) 



(67) 



In this case, it is clear that the IGE presents a power law decay, where the power is related to the cardinality I of the 
microscopic degrees of freedom characterized by correlated pieces of macroscopic information. Furthermore, the IGE 
reaches a saturation value quantified by the set {rk} (the correlation strengths). It appears that macro-correlations 
lead to the emergence of an asymptotic information geometric compression of the explored statistical macrostates 
on the configuration manifold of the model in its evolution between initial and final macrostates Q(init'a,i) ^q-j 
9 (final) (-^-j^ respectively. 



C. Ensemble of Random Frequency Macroscopic Inverted Harmonic Oscillators 



For our third example, we employ ED in conjunction with "Newtonian Entropic Dynamics" (NED) [s^]. In NED, 
we explore the possibility of using well established principles of inference to derive Newtonian dynamics from relevant 
prior information codified into an appropriate statistical manifold. The basic assumption is that there is an irreducible 
uncertainty in the location of particles so that the position of a particle is defined by a probability distribution. The 
corresponding configuration space is a statistical manifold Ms the geometry of which is defined by the Fisher- Rao 
information metric. The trajectory follows from a principle of inference, namely the method of Maximum relative 
Entropy. There is no need for additional "physical" postulates such as an action principle or equation of motion, nor 
for the concept of mass, momentum or phase space, not even the notion of time. The resulting "entropic" dynamics 
reproduces Newton's mechanics for any number of particles interacting among themselves and with external fields. 
Both the mass of the particles and their interactions are explained as a consequence of the underlying statistical 
manifold. 
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In our special application, we consider a manifold with a line element ds^ = gab (6) dd'^M^ {a,b = 1, . . . ,1) given 



ds" = [1 - $ (e)] 5ab (e) drdd\ $ (e) ^ ^^u,. {dk) 



(68) 



fc=i 



where 



(69) 



The geodesic equations for the macrovariables (s) are strongly nonlinear and their integration is not trivial. 
However, upon a suitable change of the affine parameter s used in the geodesic equations, we may simplify the 

differential equations for the macroscopic variables parametrizing points on the manifold Jvl^ =^ -^iho "^ith metric 
tensor gab- Recalling that the notion of chaos is observer-dependent and upon changing the afhne parameter from s to 
T in such a way that ds^ = 2(1 - $) dr^, we obtain new geodesic equations describing a set of macroscopic inverted 
harmonic oscillators (IHOs). In this example, the IGE S .ai) (r; wi, . . . , w;) reads 



where 



dcf 



w/) = lim lnC^(o (t; wi, . . . , w;) , 



T 

,uJi) ^ - AV i 11 (r';wi,...,Wi)rfT', 

T J ^^'iHO 



(70) 



(71) 



and 



Substituting ([7T|) and ([72]) in ([70]) . we obtain the general expression for Sj^(i) (r; wi, . . . , w;), 



(72) 



Sj^(i_l, (t; cJi , . . . , w/) =^ jim^ In \ ^ 



{3'} V 



To evaluate (|73l) we observe that Ay^(i) in ([7^ can be written as 



dr' 



(73) 



'mo t 22 \'i=l 



Since the Z-Newtonian equations of motions for each IHO are given by 

$--,^^. = o,v, = i,...,/, 

the asymptotic behavior of such macrovariables on manifold A^jHo given by 



(74) 



(75) 



(76) 



We therefore obtain 



Ay^(i) (T;a;i, . . . ,a;i) 



>oo 1 1 



j f .exp^. 



(77) 
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Upon averaging (j77p . we find 



,^i) 




UE, ] ■ exp (fir') 



(78) 



where SI — ^Wi- As a working hypothesis, we assume = Sj = S Vi, j = 1,...,Z. Furthermore, assume that 

n — > oo SO that the spectrum of frequencies becomes continuum and as an additional working hypothesis, assume this 
spectrum is Unearly distributed (Ohmic frequency spectrum). 



POhmic i^) 



f^cut-ott 



"cut-off 



UJ with / Pohmic (^) f^cut-off = ^f^, C G 



(79) 



Thus, we obtain 



Finally, substituting dSD]) into ^ yields 0, H 



l^„2/^f^^'^'exp(i^^]r) 



(80) 



5^(0 (t;wi, . . . ,iLj;) oc Qt, n 



I 



(81) 



Equation ()81|) displays the asymptotic, linear information geometrodynamical entropy growth of the generalized l- 
set of inverted harmonic oscillators and extends the result of Zurek-Paz to an arbitrary set of anisotropic inverted 
harmonic oscillators [54| in a classical information-geometric setting. In order to ensure the compactification of the 
parameter space of the system (and therefore Ms itself), it is possible to choose a Gaussian distributed frequency 
spectrum for the IHOs. With this choice of frequency spectrum, the folding mechanism required for true chaos is 
restored in a statistical (averaging over lj and r) sense. This example may be considered the information geometric 
analogue of the Zurek-Paz model used to investigate the implications of decoherence for quantum chaos. In their 
work, Zurek and Paz considered a chaotic system, a single unstable harmonic oscillator characterized by a potential 
V (x) = Lyapunov exponent), coupled to an external environment. In the reversible classical limit 

[591 , the von Neumann entropy of such a system increases linearly at a rate determined by the Lyapunov exponent. 



(^(chaotic) / \ 
'-^quantum V / 

with f7 playing the role of the Lyapunov exponent. 



(82) 



D. Regular and Chaotic Quantum Spin Chains 

In our fourth example, we apply the IGAC to study the entropic dynamics on curved statistical manifolds induced by 
classical probability distributions commonly used in the study of regular and chaotic quantum energy level statistics. 
In doing so, we suggest an information-geometric characterization of a special class of regular and chaotic quantum 
energy level statistics. 

Recall that the theory of quantum chaos (quantum mechanics of systems whose classical dynamics are chaotic) is 
not primarily related to few-body physics. Indeed, in real physical systems such as many-electron atoms and heavy 
nuclei, the origin of complex behavior is the very strong interaction among many particles. To deal with such systems, 
a famous statistical approach has been developed which is based upon Random Matrix Theory (RMT). The main idea 
of this approach is to neglect the detailed description of the motion and to treat these systems statistically, bearing 
in mind that the interaction among particles is so complex and strong that generic properties are expected to emerge. 
Once again, this is exactly the philosophy underlining the ED approach to complex dynamics. It is known that the 
asymptotic behavior of computational costs and entanglement entropies of integrable and chaotic Ising spin chains are 
very different [6^. Prosen considered the question of time efficiency in implementing an up-to-date version of the t- 
DMRG for a family of Ising spin i chains in an arbitrarily oriented magnetic field which undergoes a transition from an 
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integrable (transverse Ising) to nonintegrable chaotic regime as the magnetic field is varied. An integrable (regular) 
Ising chain in a general homogeneous, transverse magnetic field is defined through the Hamiltonian "Hicguiar (Oj 2), 
where 

n-2 n-1 

n {K,hy) = + {h-a^ + /iVJ) . (83) 

j=o j=0 

In this case, the computational cost presents a polynomial growth in time, £)('^°suiat) ^^-^ -j-, while the entanglement 
entropy is characterized by logarithmic growth, 

^regular (0, 2) = 4on Lumann ^^"^ C In T + c'. (84) 

The constant c depends exclusively on the value of the fixed transverse magnetic field intensity B± , while c' depends 
on B± and on the choice of the initial local operators of finite index used to calculate the operator space entanglement 
entropy. In contrast, a quantum chaotic Ising chain in a general homogeneous tilted magnetic field is defined through 
the Hamiltonian "Hchaotic (1, 1), where Ji is given in (j83p . In this case, the computational cost presents an exponential 
growth in time, ]j^'^^'^°^"^'> (f^ q<- g^p (/C^r), while the entanglement entropy is characterized by linear growth, 

^chaotic (1, 1) = 5^on Neumann °^ ^q'^- (^5) 

The quantity ICq is a constant, is asymptotically independent of the number of indexes of the initial local operators 
used to calculate the operator space entropy, depends only on the Hamiltonian evolution and not on the details of 
the initial state observable or error measures, and can be interpreted as a kind of quantum dynamical entropy. 

It is well known the quantum description of chaos is characterized by a radical change in the statistics of quantum 
energy levels [6l| . The transition to chaos in the classical limit of quantum systems is associated with a drastic change 
in the statistics of the nearest-neighbor spacings of quantum energy levels. In the regular regime, the distribution 
agrees with Poisson statistics while in the chaotic regime the Wigner-Dyson distribution works very well. Uncorrelated 
energy levels are characteristic of quantum systems corresponding to a classically regular motion while a level repulsion 
(a suppression of small energy level spacing) is typical for systems which are classically chaotic. A standard quantum 
example is provided by the study of energy level statistics of an Hydrogen atom in a strong magnetic field. It is 
known that level spacing distribution (LSD) is a standard indicator of quantum chaos [12] ■ It displays characteristic 
level repulsion for strongly nonintegrable quantum systems, whereas for integrable systems there is no repulsion due 
to the existence of conservation laws and quantum numbers. In [gOl, the authors calculate the LSD of the spectra of 
"Hroguiar (0, 2) and 'Hchaotic (1, !)• They find that for "Hrcguiar (0, 2), the nearest neighbor LSD is described by a Poisson 
distribution. For "Hchaotic (1, 1), they find the nearest neighbor LSD is described by a Wigner-Dyson distribution. 
Therefore, they conclude that T^reguiar (0,2) and 'Hchaotic (1, 1) indeed represent generic regular and quantum chaotic 
systems, respectively. We encode the relevant information about the spin-chain in a suitable composite-probability 
distribution taking account of the quantum spin chain and the configuration of the external magnetic field in which 
they are immersed. 

In the ME method, the selection of relevant variables is made on the basis of intuition guided by experiment; it is 
essentially a matter of trial and error. The variables should include those that can be controlled or experimentally 
observed, but there are cases where others must also be considered. Our objective here is to choose the relevant 
microvariables of the system and select the relevant information concerning each one of them. In the integrable case, 
the Hamiltonian ^regular (0, 2) describes an antiferromagnetic Ising chain immersed in a transverse, homogeneous 
magnetic field i?transversc = B± where the level spacing distribution of its spectrum is given by the Poisson 
distribution 

PPoisson (a; A I Ai^ ) = — cxp ( - — ) • (86) 
Ma V Ma/ 

The microvariable xa represents the spacing of the energy levels while the macrovariable /z^ is the average spacing. The 
chain is immersed in a transverse magnetic field which has just one component B± in the Hamiltonian T^roguiar (0, 2). 
Observe that the Exponential distribution is identified by information theory as the maximum entropy distribution 
if only one piece of information (the expectation value) is known. Thus, we translate this piece of information in the 
IGAC formalism by coupling the probability (j86p to an exponential bath p^'^P°'i<="*'ai) (^sj/i^) given by 

^(exponential) ^^^|^^ ^ ±_ / ^\ ^g^^ 

Mb V Ms/ 
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^xb\- 


exp 










^^BJ . 



where the microvariable a;^ is the intensity of the magnetic field and the macrovariable is the average intensity. 
More correctly, xb should be the energy arising from the interaction of the transverse magnetic field with the spin 

i particle magnetic moment, xb = —fj- ■ B = \—^Bcos.Lp\ where Lp is the tilt angle. For the sake of simplicity, let 

us set = 1. Then in the transverse case — Q and therefore xb = B = B±. This is our best guess and we 
justify it by noticing that the magnetic field intensity is indeed a relevant quantity in this experiment (see equation 
(|84p ). Its components are varied during the transition from integrable to chaotic regimes. In the regular regime, we 
say the magnetic field intensity is set to a well-defined value (xb) = fJ-B- Finally, the chosen composite probability 

distribution p(integrablo) 

, xbImajMb) encoding relevant information about the system is given by 

p(intcgrablc) (^^ ^ , ) = 

Again, we point out that our probability (1881) is our best guess and, of course, must be consistent with numerical 
simulations and experimental data in order to have some merit. We point out that equation (|88p is not fully justified 
from a theoretical point of view, a situation that occurs due to the lack of a systematic way to select the relevant 
microvariables of the system (and to choose the appropriate information about such microvariables) . Let us denote 
^mtograbio ^j^^ two-dimcusional curved statistical manifold underlying our information geometrodynamics. The line 
element ds^ntegrabic on Mf'^'''^^'' is given by 

'^■^hitcgrablc = '^^Poisson + '^'^Exponential = ~^dfl\ H 2~'^Ms- (^9) 

Ma Ms 

Applying the IGAC to the line element in (|89p leads to conclude polynomial growth in (2^^°^'^^^^° and logarithmic IGE 
growth [3, 111 , 

^Megrable) ™ ^^p(^,^)^c:o ^ ^Mcgrablc) .^oo ^^^^^^ ^ ^^^^ 

The quantity cig is a constant proportional to the number of Exponential probability distributions in the compos- 
ite distribution used to calculate the IGE; c[q is a constant that depends on the values assumed by the statistical 
macrovariables and fi^. Equations (|90l) may be interpreted as the information-geometric analogue of the computa- 
tional complexity £)('°suiar) entanglement entropy ^regular (0, 2) defined in standard quantum information 
theory, respectively. We cannot state they are the same since we are not fully justifying, from a theoretical standpoint, 
our choice of the composite probability ([88)) . 

In the chaotic case, the Hamiltonian Hchaotic (1; 1) describes an antiferromagnetic Ising chain immersed in a tilted, 
homogeneous magnetic field Stilted = B± B± + By , with the level spacing distribution of its spectrum given by 

the PoisSOn distribution PWigner-Dyson (2;aIMa)i 



PWigncr-Dyson (a;A|MA) 7ri2 ^ ' (^l) 



where the microvariable represents the spacing of the energy levels and the macrovariable is the average spacing. 
The chain is immersed in the tilted magnetic vector field which has two components B±^ and By in the Hamiltonian 
^chaotic (1 7 !)• The Gaussian distribution is identified by information theory as the maximum entropy distribution 
if only the expectation value and the variance are known. We translate this information in the IGAC formalism by 
coupling the probability (|9l]) to a Gaussian p^''"*'*''^'^'' {x'g\^'g^a'g) 

(Gaussian) i I \ I I \ 1 / (-^S ~ Ms) 1 

Pb {xB\fiB,aB) = ^==cxp{ , (92) 

V27rcr'J \^ 2fT^ J 

where the microvariable x'g is the intensity of the magnetic field, the macrovariable fi'^ is the average intensity of the 
magnetic energy arising from the interaction of the tilted magnetic field with the spin ^ particle magnetic moment 
and a'g is its covariance: during the transition from the integrable to the chaotic regime, the magnetic field intensity is 
being varied (experimentally). It is being tilted and its two components {B± and _By) are being varied as well. Our best 
guess based on the experimental mechanism that drives the transitions between the two regimes is that the magnetic 
field intensity (actually the microvariable iiB cos ip) is Gaussian-distributed (two macrovariables) during this change. 
In the chaotic regime, we say the magnetic field intensity is set to a well-defined value (xg) = /i'^ with covariance 
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{x'g — {x'g))'^\. Thus, the chosen composite probabihty distribution P^' 



chaotic) 



relevant information about the system is given by 



(a;^,a;slAi!4'A*s>^s) encoding 



TT (27 



2fl 



-Xj^ exp 



^ , {x'b - Mb) 



(93) 



Let us denote by three-dimensional curved statistical manifold underlying our ED model. The corre- 

sponding line element ds^Yiaotic _y\^^<=haotic) gjygjj 



^„2 _ I 2 

""'chaotic ~ ""'Wignor-Dyson 



ds 



Gaussian 



4 

17'- 



,2 '^I^'a 



-^dn 



-p^da'l. 



(94) 



Applying the IGAC machinery to the line element in ((94|). we obtain exponential growth for v^''"*"^ and linear IGE 



growth [3, [1] , 



^(chaotic) / \ T— >-oo ^ / \ o(chaotic) / \ r— foo 

^M, v) « CiG exp (/Cigt) , 5)^ ' (r) cx /Ciqt 



(95) 



The constant Ciq encodes information about the initial conditions of the statistical macrovariables parametrizing 



elements of M 



(chaotic) 



The constant A^iq, given by 



IC 



IG 



dSM, (t) 
dT 



lim 



T 



Jms (0) 



= A 



Ms 



(96) 



is the model parameter of the chaotic system and depends on the temporal evolution of the statistical macrovariables. 
It plays the role of the standard Lyapunov exponent of a trajectory and is, in principle, an experimentally observable 
quantity. The quantity J^Ms (''") the Jacobi field intensity and Xms be considered the information-geometric 
analogue of the leading Lyapunov exponent in conventional Hamiltonian systems. Given an explicit expression of /Cig 
in terms of the observables /i^, /i'^ and cr^, a clear understanding of the relation between the IGE (or /Cig) and the 
entanglement entropy (or )Cq) becomes the key point that deserves further study. Equations (195^ are the information- 
geometric analogue of the computational complexity (^-p) and the entanglement entropy iSchaotic (1, 1) defined 
in standard quantum information theory, respectively. This result requires a deeper analysis in order to be fully 
understood. One of the major limitations of our findings is the lack of a detailed account of the comparison of theory 
with experiment. This point will be among our primary concerns in future works. Some considerations may however, 
be carried out at the present stage. The experimental observables in our theoretical models are the statistical 
macrovariables characterizing the composite probability distributions. In the integrable case, where the coupling 
between a Poisson and an Exponential distribution is considered, fij^ and fig are the experimental observables. In the 
chaotic case, where the coupling between a Wigner-Dyson and a Gaussian distribution is considered, /x^, fig and a'g 
play the role of the experimental observables. We think one way to test our theory may be to determine a numerical 
estimate of the leading Lyapunov exponent Amax or the Lyapunov spectrum for the Hamiltonian systems under 
investigation directly from experimental data (measurement of a time series) and compare it to our theoretical estimate 
for Xms [s^t - However, we are aware that it may be rather difficult to evaluate Lyapunov exponents numerically. 
Otherwise, knowing that the mean values of the positive Lyapunov exponents are related to the Kolmogorov-Sinai 
(KS) dynamical entropy, we suggest to measure the KS entropy IC directly from a time signal associated with a 
suitable combination of our experimental observables and compare it to our indirect theoretical estimate for /Ciq from 
the asymptotic behaviors of our statistical macrovariables . We are aware that the basis of our discussion is rather 
qualitative. However, we hope that with additional study, especially in clarifying the relation between the IGE and 
the entanglement entropy, our theoretical information geometric characterization will find experimental support in 
the future. For these reasons, the statement that our findings may be relevant to experiments verifying the existence 
of chaoticity and related dynamical properties on a macroscopic level in energy level statistics in chaotic and regular 
quantum spin chains is purely a conjecture at this stage. 



E. Quantum Entangled Wave Packets 

As our final example, we apply the IGAC to characterize the quantum entanglement produced by a head-on 
collision between two identical (but distinguishable) spinless, structureless, non-relativistic particles of mass m, each 
represented by minimum uncertainty Gaussian wave-packets interacting via a scattering process [55| . Before colliding. 
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the two particles are in the form of disentangled Gaussian wave packets, each characterized by a width ctq in momentum 
space. The initial distance between the two particles is Rq and their average initial momenta - setting the Planck 
constant h equal to one - are T^o, respectively. After some straightforward algebra [sHi it can be shown that the 
initial (pre-collisional) two-particle square wave amplitude in momentum space is given by 



1 



■ exp 



(fcl - fco)' + (fc2 + fco)' 



(97) 



where ah 



±fco 



(^i/2>o 



{Pl/2) 



with (pi)o = Po, (P2)o = ~Po, fo is defined as in ([13]) and we 



have made use of the center of mass and relative coordinates whose conjugate momenta are defined as K = ki + k2 €z 
(—00, +00) and k = ^ (ki — k2) G (—00, +00) with — ^"Tp G (~oo, +00). 

Similarly, following [55|, and after some tedious algebra, one finds that the final (long time limit, post-collisional) 
two-particle square wave amplitude in momentum space is given by 



d(QM) 

post 



exp 



(A:i,fc2|fco,crfc„;rQM) 



'2(i--aM) 



(fcl 



2r-QM 



(fci-fco)(fc2+fco) {k2 + kaY 



'QM 



with 



S{2kl+al)R,a,<. 1, 



(98) 



(99) 



where the parameter as has dimension of length and is defined as the s-wave scattering length ^5]. It is evident 



from (1991) that tqm is non-zero and positive. Thus, / (fco) 



~ fco 



or 



fcoi<Sl 



-flg, where 



Oo = 9 (ko) ~ — '^^^^i^ ~ —koQs denotes the s-wave scattering phase shift, /i is the reduced mass /i = m/2, / (fc)is 
the scattering amplitude and L is the range of the scattering potential V given by 



V{x) = 



V, 0<x<L 
0, x> L 



(100) 



where V denotes the height (for ^ > 0; repulsive potential) or depth (for V < 0; attractive potential) of the potential. 
The quantity 9 (k) is the s-wave scattering phase shift considered around k = ko (i.e. assuming our wave-packet is 
well-localized around k — ko) and in the limit of low-energy scattering, i.e. 9 (k) <^ 1. 

We conjecture that the quantum entanglement produced by a head-on collision between two Gaussian wave packets 
are macroscopic manifestations emerging from specific underlying microscopic statistical structures. Specifically, we 
propose that Pp?J^^ (^1,^21^0,0"^^) can be interpreted as a limiting case (initial time limit) arising from a Gaussian 



probability distribution Pp 



(IG) 



dof 



Pplc^ (a;,y|Aia:>Ma,crj = ^^^^ 



exp 



2a2 



Upon setting a; -> fci, ?/ -> fc2, (x) = /i^ ^ ^^.^ = +ko, (y) = fiy 



2a2 

= —ko and a 



5(IG) 
pre 



{x,y\fi.^,fiy;a) 



P. 



pre' (fci,fc2|Mfei,Aifc2.f^) 



Pp,c^ {ki , ^2 iMfci , Aifc2 , 0-) =^ T^— J exp 



27ra-2 



2a2 



2f72 



enabling the identification 



pre 



(101) 

(Tfe^, we obtain 
(102) 

(103) 



The variances ax and ay in the random variables x and y respectively, are given by the standard definition (|13p . We 
remark that in general, a^ ay. In the present example however, it is sufficient to consider ax = ay — a. 
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We propose that -Piost^'' (^ii ^21^0, Cfc^; tqm) can be viewed as a limiting case (final or long time limit) arising from 



a Gaussian probability distribution Pp^^ (a:, yllJ-x, fJ-y', c, ?') , 

exp 



Pp^t {x,y\^l^,^ly]a,r) 



— 2 —2 _2 



27rCT2Vr 



(104) 



where the micro-correlation coefficient r is defined as in (j24l) . In the present example, the micro-correlation coefficient 
r is considered to have compact support over the line segment [0, 1), that is r G [0, 1). Upon setting /ij,^ — > +fco, 
/i^^ — > — fco and cr — >• (Tfe^ we obtain 



-Ppo?t (fci,fc2|fco,crfco;r) 



exp 



73y 



— — 



Oj. (fcl-fc°Kfc2-Hfco) 



(fc2-|-fco)^ 

— :;3 — 



27rcr2 Vl -r2 



(105) 



In this case, when both the weak correlation (r ^ 1) and the weak scattering conditions (16* (fco)| ^ 1) are satisfied, 
we obtain an excellent overlapping between ((98)) and (|105p . so that 



^p?sf^ (fci,fc2|fco,crfc„;rQM) ^ Ppo?t {k,ko,ak„;r) for r < 1, tqm < 1 and |6'(fco)| < 1, 



(106) 



assuming that fco, Ufc^, r and tqm are fixed numerical constants, and letting assume values in the neighborhood 
of ko- 

At this stage our conjecture is only mathematically sustained by the formal identities (|103p and (|106p . To render 
our conjecture physically relevant, recall that s-wave scattering can also be described in terms of a scattering potential 
V{x) and the scattering phase shift 9 (k). Integrating the radial part of the Schrodinger equation with this potential 
for the scattered wave and imposing the matching condition at a; = i for its solution and its first derivative leads to 

M 



kin cot (hnL) = fcout COt (fcouti + 0) 



with 



< X < L, 



X > L, 



(107) 

(108) 
(109) 



The quantities fi and £ are the reduced mass and kinetic energy of the two-particle system in the relative coordinates, 
respectively; km and fcout represent the conjugate-coordinate wave vectors inside and outside the potential region, 
respectively. Equation (|107l) indicates that the scattering potential V{x) shifts the phase of the scattered wave at 
points beyond the scattering region. 

Our information geometric modeling may be briefly described in the following way. The pre-coUisional scenario 
is characterized by the information geometric dynamics on the curved statistical manifold tVI^""'^"" -' of uncorrelated 
Gaussian probability distributions PpJ?^ {ki,k2\^i}^^, ^i^^^a) given in (|102p . The geodesic trajectories on for 
the non-correlated Gaussian system is given by 



(pib(r)) = Ail (r;0) = -v^2rp^tanh(Aor), 
(P2b(T)) = M2 0) = ^/pI + 2cr2 tanh {Aot) , 



(110) 

(111) 

(112) 



The post-colhsional scenario is characterized by the information geometric dynamics on the curved statistical man- 
ifold M 



(corr.) 



of correlated Gaussian probability distributions Ppost (^i; ^2|^ot cr/co5 ^) given in (jlQ5p . The geodesic 



trajectories on M 



(corr/ 



for the non-correlated Gaussian system is given by 

(pia(T)) = /ii (r; r) = - v/(l - r) {pi + 2al) tanh {A^t) , 
(p2a(T)) = /i2(T; r) = V(l - r) {pl + 2al) tanh {A^t) , 

(CTa(r)) = cr(r;r) = 



1 



2^° cosh (At)' 



(113) 
(114) 

(115) 
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where the subscript "o" denotes the initial state, the subscripts "i" and "2" denote particle 1 and particle 2, respec- 
tively; subscripts "b" and "a" denote 'before' and 'after' collision, respectively and 

Ao = — sinh-if^;^) (116) 



To V CTq / 2 \Po J 8 VPo 



, 6 

Po 



The two sets of geodesic curves comprised of {(pib(''')) , {p2h{T)) , {<^h{T))} (for the non-correlated model) and 
{{PiaiT)) I {P2a{T)) , {(Ta.{T))} (for the Correlated model) are joined at the junction r = 0: r < (before colUsion) 
for the non-correlated model and r > (after collision) for the correlated model. We recognize that the momenta 
{pvo{t)) and {pi^{t)) asymptotically converge to \/p1 + 2cr^ and — ■\/(l — r) (p^ 4- 2(t^) toward r = —00 and r = -l-oo, 
respectively (the same is true for — {p2h{T)) and — {p2a.{'T))) while {(T\,(t)) and (f7a(T)) are identical and vanishingly 
small toward r = ±cxd. Furthermore, we observe that there is continuity between (pi/2b(''')) ^nd {pi/2a{'i')') and 
between {uh{T)) and (aa^r)) at the junction, r = 

A question that now arises is how to determine the scattering phase shift in view of the fact that our statistical 
model is correlated after collision. Initially, we need to examine how correlations affect the momentum geodesic curve 
{pi/2{t))- For this purpose we define the momentum-difference curve {p{t)) = ^ [(p2(t)) — (pi(t))]. Comparison of 
the following two equations 

1 



(p(r;0)) = -[(p2b(r))-(pib(T))]- Vpl + 2^tanh(AoT), (117) 
1 



(p(T;r)) EE -[(p2a(T))-(pia(r))] = V(1^0M+2^tanh(AoT), (118) 

which follow from (jllOp . (|llip . (|113p and (I114p . indicates that at any arbitrary time r > 

(p(t;0)) > (p(T;r)), (119) 

while both (|117p and (jll8l) share the functional argument AoT. Condition (|119p implies that the correlation causes 
a reduction in the momentum for any r > (relative to the non-correlated case). This situation is analogous to 
the change in momentum caused by a repulsive scattering potential (see (|108|) and (|109p ). It is then reasonable to 
assume there exists some connection between the scattering potential and the correlation. Provided this connection 
is established, one should be able to determine the scattering phase shift in terms of the correlation via equations 
(|107p . (|108p and (|109p . In this way, one can ultimately establish a connection between quantum entanglement and 
the statistical micro-correlation. 

Recall that before collision (at the afflne time — Tq) particles 1 and 2 are separated by a linear distance Rq. Each 
particle has momenta po and —po, respectively and the same momentum spread CTq- Then from (|110p . (|llip and (jll2p 
we have 



Po = {Plh (-To)) = - (P2b (-To)) = ^P? + Sct^ tanh (AqTo) , (120) 

For arbitrary r > after collision, the system of particles 1 and 2, which initially carried momenta p^ and —Po, 
respectively at r = — Tq before collision, now carries the relative conjugate-momentum {p (r; r)) given by (jll8l) due to 
the correlation. With nonvanishing micro-correlation the wave-packets experience the effect of a repulsive potential; 
the magnitude of the wave vectors (or momenta) decreases relative to the corresponding non-correlated value. One 
may rewrite (fTOTl) . (fTOHl) and ([T09| as 

kr cot {kj.L) = fco cot {koL + 9o) , (122) 

with 



fc. = ^^^fHZ), o<x<L, (123) 

h 

fco = x>L, (124) 
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where kr and fco represent the wave vectors with and without the correlation, respectively. The connection between 
the correlation and the scattering potential can be established by combining (|123p and (11241) . From (|119p one finds 
that the correlation renders 



ko — > kr = VI - rko. (125) 
Then using ()123|) . ()124p and ()125|) . we determine the scattering potential, 

V^r£^r^=r^. (126) 

Equation (jl26p clearly establishes a connection between the correlation coefficient and the scattering potential: the 
correlation coefficient is the ratio of the scattering potential to the initial relative kinetic energy of the system. From 
(jl26p it is evident that our interaction potential is repulsive, i.e. y > since we consider non-negative micro- 
correlations, r e [0, 1). 

With the potential determined, one can determine the scattering phase shift 9o (for low energy s-wave scattering) 
by combining equations and the result being, 

une^'^^'^'-'^'e^^-!^. (127) 

By means of p26p and (I127|) we can express the scattering phase shift in terms of the scattering potential 

2^iVkoL^ 2fiVpoL^ 



(128) 



which is in agreement with j67| . This is the first significant finding that allows to state that our conjecture is also 
physically motivated. 

As the scattering potential has been determined, so too can the scattering amplitude be determined. To this end, 
we write 

e*^°sin6'o _ 6'o _ .,„„^ 
f(ko) = 7 ~ T ~ -'^s (129) 

for low energy s-wave scattering, koL — p^L/h <C 1. Thus, we finally obtain the scattering cross section: 

E = 47r|/(A:o)| ~ = ^4^3. (130) 

In order to properly analyze entanglement, the entanglement entropy obtained from the long time limit post-collisional 
wave function is required. In most cases however, this must be performed numerically. Thus, to approach the problem 
analytically and simultaneously gain insights into the problem, it is convenient to make use of the linearized version 
of the entropy of the system, i.e. of the purity of the system [55|- The purity function is defined as 

7'-Tr(p^), (131) 

where = Tr^ {pab) the reduced density matrix of particle A and is the two-particle density matrix 
associated with the post-collisional two-particle wave function. For pure two-particle states, the smaller the value of 
V the higher the entanglement. That is, the loss of purity provides an indicator of the degree of entanglement. Hence, 
a disentangled product state corresponds to V = 1. Under the assumption that the two particles are well separated 
both initially (before collision) and finally (after collision), and further assuming that the colliding Gaussian wave 
packets are very narrow in momentum space ((7^^ <C 1 such that the phase shift can be treated as a constant 6 (fco)), 
it follows that the purity of the post-collisional two-particle wave function is approximately given by [55|] 

7^ = 1 - 8 (2fc2 + cr2J ^^as -I- O {al) . (132) 

Employing the scattering cross section E = 47r |/ (fco)|^ ~ 47rr- k^L _ istt^^^Y ^ w ^iral, we may express the purity 
in an alternative manner, namely, 

4 (2fc2 + al ) i?oVE 



23 



Equations (|132[) and (|133l) above demonstrate how the entanglement can be measured from the loss of purity by use 
of the scattering length or cross section. By combining (|132p and the square of (|129p we find the purity 



1 - 



8rfc2 (2fc2 



RoL^ ^ mfiV (2fc2 



RoL^ 



(134) 



Equation ()134[) implies that the purity V can be expressed in terms of physical quantities such as the scattering 
potential height V and the initial quantities fco, fo and Ro via (|130p and (jl26l) . This is the second significant finding 
obtained within our hybrid approach (quantum dynamical results combined with information geometric modeling 
techniques) that allows to explain how the interaction potential height V and the incident particle energies^ control the 
strength of the entanglement. The role played by r in the quantities V and V suggests that the physical information 
about quantum scattering - and therefore about quantum entanglement - is encoded in the statistical correlation 

coefficient, specifically in the covariance term Cov(fci,A:2) =^ {kik2) — (fci) (A;2) appearing in the definition of r. 

The correlation coefficient r can now be expressed in terms of the physical quantities such as the scattering potential, 
the scattering cross section and the purity. Solving equations (I126p . (|130p and (|134p for r, we obtain 



2nV 2fiV 



3ys 



2v^fc2L3' 

3(1-7^) 

In view of (IM)) . (|130l) and (I136p . one obtains the following relation: 

V ^ 4h^kt {2kl + aj^) R, 
3^ 



(135) 
(136) 
(137) 

(138) 



which indicates the uniform scattering potential density is solely determined by the initial conditions of the given 
system. 

From (|117p . (I118p and (|119p it is observed that for the micro-correlated Gaussian system considered here, more time 
is required to attain the same momentum value compared with the non-correlated Gaussian system. For example, 
in order to attain the same value as the initial momentum po, the non-correlated system and the micro-correlated 
system would require time intervals Tq and r* , respectively, where 



Po 
Po 



Vp^ 



2a-2tanh (AqTo) , 



r) (j32 + 2cr2)tanh(AoT*). 



Combining (IT^ and ([Till)) , we obtain 

tanh (AoT^) = (1 
Rewriting and expanding both sides of (I14ip . we have 

l-2e-2^°"* +0(e-4^°"*) = (1- 
Rounding (|142p off and arranging terms we obtain 



-1/2 



tanh (AqTo) 



-1/2 



[1 - 2e 



O e 



)] 



(1 



-1/2 



1 r 
2 



(1 



-1/2 



The first term on the right hand side of (|143p can be approximated to 1 since 
However, r in the second term should not be disregarded in the same manner because (1 — r 
[ir -I- O (f^)] e^^°'^° is not negligible. Therefore, we may rewrite (|143p as 



,-1/2 



(139) 
(140) 



(141) 
(142) 

(143) 
r^) and r ^ 1. 



1 



o2Ao 



-2AoA 



1 - 



(1-r) 



-1/2 



- 1 



(144) 
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where A = — Tq represents a new quantity that we term " prolongation" , and 



exp 



O 



(145) 



for ^ ^ 1 due to (jll6p . The quantities and Tq are the temporal intervals required for a particle to reach the 
same value of momentum fco from in the post-coUisional scenario, in presence and in the absence of correlations r, 
respectively. From ()144|) we find 



A(A:o,ao,r) oc ln{l - (1 - r) 



-1/2 



(146) 



Here, we can find the upper bound value of r by means of (|146p and (|145p . 



r < 



'7a 



(147) 



The prolongation serves to quantify the time required by a micro-correlated system - relative to a corresponding 
non-correlated one - to attain the same momentum value (relative to the same initial reference time) . The occurrence 
of a non- vanishing prolongation is in fact due to the existence of micro-correlations and therefore, due to the existence 
of scattering phase shifts. In other words, in the absence of scattering there is no time difference. This can be stated 
in yet another way as follows: "The prolongation encodes information about how long it would take an entangled 
system to overcome the momentum gap (relative to a corresponding non-entangled system) generated by the scattering 
phase shift. The entangled system only attains the full value of momentum (i.e. the momentum value as seen in 
the corresponding non-entangled system) when the scattering phase shift vanishes. For this reason, the prolongation 
represents the temporal duration over which the entanglement is active" . We observe that the entanglement duration 
can be controlled via the initial parameters fco, ctq and the correlations r (therefore via the incident particle energies 
and the scattering potential due to (ll35p V Moreover, we notice that in the absence of correlations (i.e. r — )■ 0), 
A — >■ 0. It is anticipated that the maximum duration would be obtained when r is the greatest and the ratio (To/fco 
is the smallest. 

The line element ds^ = gab (6) M'^dd'' (a, 6 = 1,2, 3) on A^i™"'^ is given by 



ds 



1 — r' 



2r 



1 



jdfi^dfiy + Ada' 



(148) 



where we consider positive micro-correlation coefficients r G (0, 1). In this geometry, the components of the sectional 
curvature are given by 



ICu 



(149) 



The Ricci scalar 72.^(corr.) and Weyl Projective Wabcd curvatures are given by 



7^^(corr.) = -- = 7^^(u„co„.) aud Wabcd = 0, 



(150) 



respectively. The fact that Wabcd — implies the manifold AA^s°^^'^ is isotropic. It is known that the anisotropy of the 
manifold underlying system dynamics plays a crucial role in the mechanism of instability. In particular, fiuctuating 
sectional curvatures require also that the manifold be anisotropic. 

It can be shown j2^] that the JLC-equation ([M)) on TWi™'^' -' reduces to 



QJ, 



= 0. 



(151) 



where Q = ^ -A^ < and 

assumes the form 



Ml 



is defined in (|37)) . Since Q < 0, unstable solutions of equation (|15ip 



(r) 



(152) 
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(eorr.)(^) 

where u (0) = = 



Recalling the definition of the hyperbolic sine function sinhx = i (e^ — e it 



is clear that the geodesic deviation on M 



(corr.) 



field intensity 



is described by means of an exponentially divergent Jacobi vector 

2 



Mi""" '' ' 



a classical feature of chaos. In order to evaluate (1551) we use (|152p to find 



'(0) 



sinh 



^g- oiiiii -Qr) and 
exponents A^(cori.) = liniT 



w2 (0) cosh^ 



Thus, for the case being considered, the Lyapunov 



\{l-Q) e^v^- 



= 2\/—Q. Therefore, it follows that 
2^/^ = 2A„ > 0. 



(153) 



From ()153p we observe the following points: the classical chaoticity does not depend on the statistical correlation, 
i.e. A^(corr.) = (non-con.) = \m, = 2 , aud the Lyapunov exponents can be determined solely from the initial 
conditions (see equation (IllGI) ). 

Yet another finding uncovers an interesting quantitative connection between quantum entanglement quantified by 
the purity V in (jl34p and the information geometric complexity (IGC) of motion on the uncorrelated and correlated 

curved statistical manifolds ) q^yA jWs°"'^ respectively. The information geometric complexity as defined in 

[68] represents the volume of the effective parametric space explored by the system in its evolution between the chosen 
initial and final macrostates. The volume itself is in general given in terms of a multidimensional integral over the 
geodesic paths connecting the initial and final macrostates. For additional details see [6^. Here, omitting technical 
details and following the works presented in [24, 69, 73|, one finds 
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where Cjq°" -* denotes the IGC on M.^s°"''' ■ 



Similarly, cj^"'^™' ' 



^|™" ) j-^ _^ Q-j represents the IGC on M. 



(154) 



(uncorr.) 



As a side remark, we point out that (jl54p confirms that an increase in the correlational structure of the dynamical 
equations for the statistical variables labelling a macrostate of a system implies a reduction in the complexity of the 
geodesic paths on the underlying curved statistical manifolds [23l.l69l|. In other words, making macroscopic predictions 
in the presence of correlations is easier than in their absence. 

The technical details that will be omitted in what follows may be found in [23] . By direct computation, the IGE is 
found to be 



S 



M 



(corr.) (T;r) ^=°° Xm,t - ln(AA^,T) + ^ In 



For non-correlated Gaussian statistical models the IGE is given by S 



(t;0) = S 
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(r; r — > 0). Observe 



that in contrast to the macro-correlated case, the IGE for the micro-correlated case presents linear growth in the 
affine temporal parameter t. Both the IGC and the IGE decrease in presence of micro-correlations. In particular, the 

IGC decreases by the factor \f^^ < 1 for r > whereas the IGE decreases by i In (t^) < for r > 0. With the 

:.) (t; r) in hand, we make the following interesting observations. From (|154p we 



quantities cj™" ' (r; r) and S 



find 
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Combining (|134p and (|154|) it follows that 



(159) 
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where the dimensionless coefScient ry^ reads 

Tjc = lkl{2kl + al)RoL\ (160) 

From (|159[) it is evident that the scattering-induced quantum entanglement and the information geometric complexity 
of motion are connected. In particular, when purity approaches to unity (entanglement- free scenario), the difference 
between the correlated and uncorrelated information geometric complexities approaches zero. 

VI. CONCLUSIONS 

In this article, we have introduced a theoretical construct that allows us to describe the macroscopic behavior 
of complex systems in terms of the underlying statistical structure of their microscopic degrees of freedom through 
statistical inference and information geometry methods. We reviewed the Maximum relative Entropy (MrE) formalism 
and the theoretical structure of the information geometrodynamical approach to chaos (IGAC) on curved statistical 
manifolds A4s- Special focus was devoted to the description of the roles played by the sectional curvature ICmsj 
the Jacobi field intensity Jms ^^.d the information geometrodynamical entropy Sms (IGE) as information geometric 
indicators of chaoticity (complexity) . Four applications of these information geometric techniques combined with ME 
methods were presented. 

First, we studied the chaotic behavior of a Gaussian statistical model describing an arbitrary system of I uncorre- 
lated degrees of freedom and found that the hyperbolicity of the non-maximally symmetric 2/-dimensional statistical 
manifold Ms underlying such a Gaussian model leads to linear IGE growth and to exponential divergence of the Jacobi 
vector field intensity [22]. Second, we studied the asymptotic behavior of the dynamical complexity of the maximum 
probability trajectories on Gaussian statistical manifolds in presence of correlation-like terms between macrovari- 
ables labeling the macrostates of the system under investigation. In presence of correlation-like terms, we observed 
a power law decay of the information geometric complexity at a rate determined by the correlation coefficient. We 
also presented an information-geometric analogue of the Zurek-Paz quantum chaos criterion of linear entropy growth. 
This analogy was motivated by studying the information geometrodynamics of an ensemble of random frequency 
macroscopic inverted harmonic oscillators. The IGAC was also employed to study the entropic dynamics on curved 
statistical manifolds induced by classical probability distributions commonly used in the study of regular and chaotic 
quantum energy level statistics. In doing so, we suggest an information-geometric characterization of regular and 
chaotic quantum energy level statistics. 

Finally, the IGAC was used to describe the scattering-induced quantum entanglement between two Gaussian wave- 
packets. The IGAC was used to analyze our specific two-variable micro-correlated Gaussian statistical model. The 
manifolds ^Ai'^°"''' and were used to model the quantum entanglement induced by head-on scattering 

(in the s-wave approximation) of two spinless, structureless, non-relativistic particles, each represented by minimum 
uncertainty wave-packets. Equation (jl28p allowed to connect the entanglement strength - quantified in terms of purity 
- to the scattering potential and incident particle energies f (|135l) and (11341) ). It was also found to be possible to relate 
the statistical entanglement duration A to the scattering potential V (x) and incident particle energies £ f (|135p and 
(I146p ). Recall that the prolongation A was defined as the time required for the observed momentum difference between 
a correlated and corresponding non-correlated system to vanish. The prolongation encodes information about how 
long it would take an entangled system to overcome the momentum gap generated by the scattering phase shift. The 
entangled system only attains the full value of momentum (i.e. the momentum value as seen in the corresponding 
non-correlated system) when the scattering phase shift vanishes. For this reason, the prolongation represents the 
temporal duration over which the entanglement is active. 

The micro-correlation coefficient r, a quantity that parameterizes the correlated microscopic degrees of freedom of 
the system, can be understood as the ratio of the potential to kinetic energy of the system. When r ^ the wave- 
packets experience the effect of a repulsive potential; the magnitude of the wave vectors (momenta) decreases relative 
to their corresponding non-correlated value. The upper bound value of r depends on po and Co in such a manner 
that r increases as Po decreases. This result constitutes a significant, explicit connection between micro-correlations 
(the correlation coefficient r) and physical observables (the macrovariable po). For r values close to its upper bound, 
the prolongation A becomes infinitely large. On the other hand, with r vanishing (i.e., no micro-correlation) A is 
identically zero. With r fixed however, the prolongation A depends on po and Uo. Thus, the prolongation A can 
be controlled by the initial conditions Po and CTq as well as r. Maximal prolongation occurs when r is greatest and 
the ratio <7o/Po is smallest. For small initial r and Po, A would be correspondingly small, suggesting that for such 
scenarios quantum entanglement is transient. Furthermore, a quantitative relation between quantum entanglement 
(purity) and the IGC (|159p was uncovered. 
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The complexity of geodesic paths on A^i'^"'^'^ -' and _A/J ' -was characterized through the asymptotic computation 

of the IGE and the Lyapunov exponents on each manifold. The Lyapunov exponents in both cases were found 
to be the same positive definite constant, A^ccon.) = A^(uncorr.) = A^i^ = ^A^ > 0. The IGE 5^(cor,.) (t; r) in 
presence of micro-correlations assumes a smaller initial value relative to the non-correlated case 5^(uncorr.) (r; 0) while 
the growth characteristics of both correlated and non-correlated IGEs were found to be similar. Specifically, the 
larger the micro-correlation (i.e. the closer r is to 1) the lower the initial value of the IGE. Thus, the stronger 
the initial micro-correlation, the larger the gap between 5^(corr.) (r — 0; r) and 5^(uncorr.) (r — 0;0). This implies 

that iS^(corr.) (T;r) < iS^(uncorr.) {t O) . Whcu micro-correlations vanish (i.e. when r = 0), we obtain the expected 
result iS^(cori.) (r; 0) = iS^(uncorr.) (t; 0). The appearance of micro-correlation terms in the elements in the Fisher-Rao 

information metric leads to the compression of C^q"'^ (t; r) by the fraction ^ and thus, to a reduction of the 
complexity of the path leading from ©(initial) ©(final) ^ 

We emphasize that at this stage of development, IGAC remains an ambitious unifying information-geometric 
theoretical construct for the study of chaotic dynamics with several unsolved problems. However, based on our 
findings, we believe it provides an interesting, innovative and potentially powerful way to study and understand the 
very important and challenging problems of classical and quantum chaos through statistical inference and information 
geometric techniques. In order to clarify the relationship between our IGE and conventional measures of complexity, 
for instance the topological entropy (the supremum of the Kolmogorov-Sinai metric entropy) , further investigation is 
required. Finally, we believe our IGE may play an important role in both classical and quantum information science, 
but at this moment this remains a conjecture 0, Q • 

One final remark: In this work, we have used information geometric techniques and inductive inference methods for 
tackling some computational problems of interest in classical and quantum physics. Specifically, we have provided an 
information geometric characterization of the complexity of dynamical systems in terms of their probabilistic descrip- 
tion on curved statistical manifolds. The maximally probable trajectories of the system arise through implementation 
of a Principle of Inference, the ME method and the constructed indicators of complexity of such trajectories are 
defined in terms of asymptotic temporal averages. This leads to the following consideration. Most standard charac- 
terizations of chaoticity are based on the Boltzmann's ergodic hypothesis which relies on a "frequency" interpretation 
of probabilities where concepts such as coarse graining (or, randomization) appear. Such artificial concepts are not 
needed in our information geometric characterization based on ME methods [l^ and they may even lead to incorrect 
results [13, [7l|. Does this mean that ergodic theory is unnecessary for our purposes [73]? Certainly, ergodic theory 
is explicitly used in the characterization of chaos in terms of the Kolmogorov-Sinai dynamical entropy and Lyapunov 
exponents [t^ . It appears that it is implicitly used in the formulation of the IGAC and the IGE by way of focusing 
on the above mentioned asymptotic temporal averages. Thus it remains to be seen whether or not the IGAC and the 
IGE are "special cases" of a broader picture. 
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